Single Dish Polarization Calibration 



Simon Johnston 



School of Physics, University of Sydney, NSW 2006, Australia 

Abstract 

Using the formalism of Hamaker et al. (1996), I derive a method for the 
polarization calibration of observations made with a single radio telescope. This 
method is particularly appropriate for observations of pulsars, where the sign and 
magnitude of the circular polarization are useful for understanding the emission 
processes at work. I apply the method to observations of PSR J1359-6038 made 
using the multibeam receiver on the Parkes radio telescope. 

1 Introduction 

The polarization properties of astronomical sources are interesting for a number of 
reasons. They can be used, for example, in determining the physical processes at work; 
thermal sources are not polarized, whereas non-thermal sources may have a significant 
degree of polarization. The rotation measure of a source can be used to determine 
magnetic field strengths; for example pulsar rotation measures are used to derive mag- 
netic maps of the Galaxy. However, the degree of polarization in astronomical sources 
is often small, particularly circular polarization, and so careful calibration is necessary 
before the results can be interpreted. 

Recent papers on polarization calibration of radio telescopes include those by Hamaker 
et al. (1996) and Britton (2000). The Hamaker et al. (1996) paper lays down the math- 
ematical foundations which are subsequently used by Sault et al. (1996) to construct 
a method for instrumental calibration of a radio interferometer. Their method is rou- 
tinely applied to data from the Australia Telescope Compact Array and incorporated 
in the software package miriad. In this paper, I will use the approach of Hamaker 
et al. (1996) to develop a method for calibration of single dish data. This method is 
particularly useful for observations of pulsars. In the final section of the paper I apply 
the method to observations of the pulsar PSR J1359-6038 using the 20-cm multibeam 
receiver on the Parkes radio telescope and discuss the results. 



2 Representation 



Following Hamaker et al. (1996) I describe the propagation of an electromagnetic wave 
in the xyz coordinate system by 
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where the axis of propagation is the z axis, and e x and e y are complex. The coherency 
properties of the electric field can be expressed in xy coordinates by the coherency 
matrix (Born & Wolf 1964) 
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where the * denotes the complex conjugate. The electromagnetic field is converted to 
electric voltage in the feed of a radio telescope which, in the case of linear feeds, consists 
of two input probes aligned along x and y. A correlator multiplies and averages these 
voltages to produce the voltage vectors 
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Generally, one is interested in the (true) Stokes parameters, /, Q, U and V which in 
combination form the Stokes vector e s . One can express e s in terms of the electric 
field vector, e, by 
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The propagation of the radiation from the astrophysical source of interest to the 
output from the correlator can be described in the following way (Hamaker et al. 
1996). Let the effects of Faraday Rotation (both through the interstellar medium and 
the ionosphere) and parallactic angle rotation be combined into the matrix R. We let 
the matrix F represent the feed response (including gains, phases and leakage terms). 
Hence the voltages at the two probes can be described by 
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and the correlator voltage vector is given by 

v = (F <g> F*) (R <g> R*) T^ 1 e s 
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where ® denotes the Kronecker (or outer) product (Hamaker et al. 1996). The rotation 
matrix, R, and the feed response matrix, F, are given by Hamaker et al. (1996) as 
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where each of B, C, G and E are complex terms. It is important to note that Hamaker 
et al. (1996) have shown that such a formalism does not involve any approximations, 
although often, as below, the matrices are then expanded only to first order. 

In this formalism, however, unlike that of Britton (2000), there is no direct physical 
association of the variables B, C, G and H. B and C can be thought of as 'leakage' 
terms; they describe the leakage of the opposite polarization into either receptor. The 
real parts of G and H relate to the sensitivity of the two probes and the difference in 
the imaginary parts of G and H relates to the phase offset between the two channels 
(often called 'instrumental phase' for a pair of linear receptors). However, this is not 
strictly correct, as G and H also mixed with the leakage terms (Hamaker et al. 1996). 

Expanding equation |6] and ignoring (for now) the effects of rotation we thus obtain 
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where the scaling factor 1/2 arises from the determinant of T. Multiplying the two 4x4 
matrices together we arrive at 
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3 Measured Stokes Parameters 



Using equation | above and defining the measured Stokes parameters in the case of 
linear feeds by 
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we can express the measured Stokes parameters in terms of the true Stokes parameters 
as follows: 
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U(gih + g 2 b 2 - hd - h 2 c 2 ) + V(gib 2 - ^2^1 + hid - /i 2 c 2 ) 
Qm = \l(gl + gl-h\-hl) + ^Q(gl + gl + h\ + hl) + 

U{g 1 b x + g 2 b 2 + hid + h 2 c 2 ) + V{g\b 2 - g 2 b x - h x c 2 + h 2 d) 
U m = + h 2 b 2 - - g 2 c 2 ) + Qi-gid- g 2 c 2 - h x bi - h 2 b 2 ) + 

U (g x hi + g 2 h 2 ) + F(#i/i 2 - #2/11) 
Kn = I(9iC2-g2Ci + h 1 b 2 -h 2 b 1 ) + Q(g 1 c 2 + g 2 d-h 1 b 2 + h 2 b 1 ) + 

U {-gih 2 + griii) + ^(^1/11 + 5-2^2) 

where all second order terms in B and C have been omitted. Somewhat unconven- 
tionally I have defined each complex term above by e.g. G = g\ + ig 2 rather than 
G = gie m to simplify the notation. Note that these parameters are assumed to be 
time-independent but will not, in general, be frequency independent. In a perfect sys- 
tem one has b\ = b 2 = 0, c% = c 2 = 0, gi = hi = 1 and g 2 = h 2 = and hence 
the measured Stokes parameters are identically equal to the true Stokes parameters 
(ignoring the effects of rotation of Stokes Q into U along the line of sight). 

If the gains and instrumental phase have been solved for by some other method, as 
is often the case, (e.g. through observations of a polarized calibration signal injected 
directly into the feed), then gi — hi — 1 and g 2 = h 2 = and equation [II] simplifies 
to: 

I m = I + U(b 1 -d) + V(b 2 + c 2 ) (12) 
Qm = Q + U(b 1 +d) + V(b 2 -c 2 ) 
U m = U + I(h-d)-Q(bi + d) 
V m = V + I(b 2 + c 2 )-Q(b 2 -c 2 ) 

We also note, as Britton (2000) has done, that a relationship holds between the 
measured and true Stokes parameters via 

%-C? m -Vl-VZ = K{I°-Cf-V*-V*) (13) 

where K is time-invariant and is related to the determinant of the matrix in Equation 
11. This expression, the so-called invariant interval, is not frequency- invariant however, 
as it depends on the (frequency dependent) terms in the F matrix. Britton's idea of 
using the invariant interval for timing purposes has been used to great effect by van 
Straten et al. (2001). 

4 Effects of feed rotation 

It is sometimes stated that performing a (short) observation, rotating the feed through 
90° and observing again and summing the two observations has the effect of cancelling 
the effect of the leakage parameters. However, this is not the case, even to first order. 
The effect of feed rotation effectively converts Q to — Q and U to — U. To sum the 
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two observations one essentially averages the sum of the two values of I m and V m and 
averages the differences of the two Q m and U m measurements. Assuming the gains are 
correctly calibrated, then using equation above, one obtains: 

I m = I + V{b 2 + c 2 ) (14) 

Q m = Q + Uih + a) 

U m = U-Q{b x + Cl ) 

V m = V + I(b 2 + c 2 ) 

and the total linear polarization is thus 

L m = Ly[\ + (h + Cl ) 2 (15) 

These last two equations show that the measured Stokes parameters are now a better 
approximation to the true Stokes parameters than a single observation would yield. 
However, even though the error in I m and L m are likely to be small, the error in V m 
can potentially be large if V is small and b 2 + c 2 is significant. The error in the position 
angle of the linear polarization will be of order tan _1 (&i + c\). 

A further potentially interesting measurement is to subtract V m from one observa- 
tion with that from an observation where the feed is rotated through 90° . In this case 
one obtains 

V m = -Q(b 2 -c 2 ) (16) 
If Q is large, one can then obtain b 2 — c 2 directly. 



5 Effects of observing a polarized calibration signal 

It is often the case that a cal probe is included in the receiver package. This probe 
usually injects a signal between the two (linear) probes at an angle of 45° to both. 
The cal thus has I = U and Q — V — 0. Imagine now that the two probes have 
identical gains and there is no instrumental rotation (i.e. we can set gi — hi — 1 and 
g 2 = h 2 = 0). Then, if the leakage terms are present, observations of the cal will yield 
a signal 

v xx * = htJ(l+2bx) (17) 

Vyy* = -U(l-2 Cl ) 

V m = U(b 2 +C 2 ) 

The first two expressions look like gain difference between the two probes, and the 
presence of V looks like instrumental rotation of U into V. Hence, even though we 
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started with a system with perfect gains and phases, we derive using equation IT2 
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and hence the system looks impure! One then re-observes the cal (or any other pure 
U signal) with these terms in equation 11 and one obtains: 
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This implies that one does not, in general, measure a pure U signal in spite of the 
calibration procedure. 



6 Circular Feeds 



The Parkes radio telescope has receivers with linear feeds. However, many other tele- 
scopes have circular feed systems. To apply this method to receivers with circular feeds 
one needs to change the Stokes vector e s to 



and equation 10 to read 
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The algebra can then be worked through to obtain expressions for the measured Stokes 
parameters in terms of the true parameters as in equation O. 



6 



7 Application 



From equation |TT] it can be seen there are 4 measured quantities (the Stokes parameters) 
and a total of 12 unknowns. These 12 unknowns comprise the 4 leakage parameters 
(&i, 62, ci, C2), the 4 gain terms (pi, #2, ^1, ^2), the total intensity, linear and circular 
polarization of the source and an unknown angle <p which deals with rotation along the 
line of sight such that Q = Lcos0 and U = Lsiiufi. This angle is made up of the intrinsic 
position angle of the source, the rotation measure through the interstellar medium and 
ionosphere, the parallactic angle and the feed angle (s) with respect to some reference 
frame. Of these angles, only the ionospheric rotation measure and the parallactic 
angle are time variable. Unless some independent measurement can be made of the 
ionosphere, this introduces unwanted noise into the calibration matrices. Ignoring 
this, rotating the feed relative to the sky (for example as naturally occurs with a feed 
mounted on an alt-az telescope) can be used to determine the 12 unknowns. Provided 
we have sufficient independent measurements of the Stokes parameters at a variety of 
parallactic angles, we can, in principle, determine all the unknowns. Sault et al. (1996) 
show that observations of an unpolarized source will yield only 9 of the unknowns, but 
that observations of a strongly linearly polarized source can determine all 12. 

Pulsars make ideal sources for (self-) calibration as they are generally highly linearly 
polarized and have a reasonable degree of circular polarization. Furthermore, the 
polarization changes across the pulse profile and one can thus use the information in 
multiple phase bins at a given pointing to help solve for the unknowns. The main 
drawback to using pulsars as calibrators is that they scintillate which implies that the 
total intensity, /, can vary significantly on timescales of order minutes. The fractional 
polarization is not affected by scintillation, however, and therefore to overcome this 
problem one is forced to normalise the measured Stokes parameters by I m at the expense 
of being unable to measure the true value of /. Unfortunately this procedure introduces 
errors as, of course, I m is not equal to /. In practice these errors will be small unless 
the source is highly polarized and bi — c\ is large. 

A previous method for calibration of pulsar signals was given in the Appendix of 
Stinebring et al. (1984) and is still in use (e.g. Weisberg et al. 1999). However, their 
method involves a number of simplifications and assumptions which, in my view, are 
incorrect in some cases and which are, in any case, no longer necessary to make. Any 
least squares fitting algorithm can easily solve for 12 unknowns with ~100 data points 
in a matter of seconds and the full equations should thus be used. 

Data on PSR J1359-6038 were collected on two separate occasions in July and 
Novmeber 2000. Each observing session lasted approximately 10 days and we obtained 
65 and 66 independent observations covering the whole parallactic angle range available 
(approximately ±94° ). These data were taken using the 64-m Parkes telescope using 
the centre beam of the multibeam receiver at a central observing frequency of 1318 
MHz and with a total bandwidth of 128 MHz. At the start of the observing session 
observations were made of the flux calibrator Hydra A to obtain the system equivalent 
flux density. Then, a 90 sec observation of the cal signal was made followed imme- 
diately by a 3 min observation of the pulsar. The pulsar observation was calibrated 
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for differential gain and phase in the two probes based on the results from the cal 
observation and flux calibrated from the observations of Hydra A. Pulse profiles were 
formed in each of I m , Q m , U m and V m with 256 phase-bins per profile and 8 frequency 
channels across the 128 MHz total bandwidth. Q m , U m and V m were then normalised 
by I m . 

Figure 1 shows the measured (normalised) Stokes parameters as a function of par- 
allatic angle from data taken in November 2000 from the peak of the pulse profile in 
one of the frequency channels centered at 1365 MHz. It can clearly be seen that the 
circular polarization, V m , varies significantly with parallactic angle and that the vari- 
ations are in phase with the variations in Q m . The deviation in V m is very large with 



a peak-to-peak amplitude of about 0.2. The implication of this (from equation [12]) is 
that b 2 — c 2 is of this order, given that Q m is about unity. 
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Parallactic Angle (degrees) 



Figure 1: Measured Stokes parameters for PSR J1359-6038 as a function of parallactic 
angle. Here, L m is denoted by circles, Q m by squares, U m by triangles and V m by 
crosses. The lines denote the best fit to the data after solving equation ITT. 



It was shown in section 5, that observations of the cal induces subsequent errors 
in the observed Stokes parameters. These errors can be removed by (re-) solving for 
the gains and phase terms (although nominally the cal observations are used to set 
these) using equation |TT[ Practically, however, I set g± = 1 and hi = and leave 
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them fixed; g 2 and h 2 then measure the relative gain and phase differences between the 
two probes. I then used a least squares fitting algorithm (in this case the Levenberg- 
Maquart Method; LMM) to minimise the residuals of (\Q m - Q\ + \U m - U\ + | V m - V\) 
and solve for the 9 unknowns from equation |11| given the data (remembering that I m 
is fixed at 1, all the observed Stokes parameters are normalised by J m and g\ = 1, 
hi = 0). The result of the fitting is shown in the Figure. At this particular frequency, 
bi = 0.0, b 2 = —0.05, Ci = 0.02 and c 2 = 0.12. Given this and the discussion in Section 
4, rotating the feed through 90° and observing again still induces a leakage of I into V 
of order 7%, about the same size as the true V signal in many pulsars. As pointed out, 
an independent measurement of 62 — C2 can be obtained in this way, and confirms the 
measurements derived above. 




-100 -50 



Parallactic Angle 

Figure 2: Rotation measure of PSR J1359-6038 as a function of parallactic angle before 
polarization calibration. The error bars are formal errors from the fitting process used 
to determine the RM. 



At this point, all the parameters of the matrix form of Equation 11 have effectively 
been determined. The true Stokes parameters can then be obtained by multiplying the 
observed Stokes parameters with the inversion of this matrix. This can be achieved 
with e.g. the numerical recipes routine lubksb. 

Figure 2 shows the rotation measure of the pulsar as a function of parallactic angle 
prior to polarization calibration. There are clear systematics in the data with sinusoidal 
variations of ~10 rad m -2 about a mean of 35 rad m~ 2 . These are due to the frequency 
dependency of the leakage terms B and C. Figure 3 shows the effects of the calibration. 
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The mean rotation measure is now 16 rad m 2 and the systematic effects are largely 
absent. 
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Figure 3: Rotation measure of PSR J1359-6038 after polarization calibration. 



8 Conclusions 

Using the excellent Hamaker et al. (1996) paper as a mathematical foundation, I 
have derived the matrices linking the true Stokes parameters to the observed Stokes 
parameters for single dish observations. It is then a computationally trivial task to 
solve for the terms in this matrix by least squares minimisation and thus provide a 
polarization calibration of the feed and receiver system. 
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